Motion in a rocked ratchet with spatially periodic friction 
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We present a detailed study of the transport and energetics of a Brow- 
nian particle moving in a periodic potential in the presence of an adiabatic 
external periodic drive. The particle is considered to move in a medium with 
periodic space dependent friction with the same periodicity as that of the po- 
tential but with a phase lag. We obtain several results, most of them arising 
due to the medium being inhomogeneous and are sensitive to the phase lag. 
When the potential is symmetric we show that efficiency of energy transduction 
can be maximised as a function of noise strength or temperature. However, in 
the case of asymmtertic potential the temperature may or may not facilitate 
the energy conversion but current reversals can be obtained as a function of 
temperature and the amplitude of the periodic drive. The reentrant behaviour 
of current can also be seen as a function of phase lag. 
PACS number(s): 05.40+j, 05.60+w, 82.20Mj 



I. INTRODUCTION 

The search for the possibihty of unidirectional motion in a periodic system without the 
apphcation of any obvious bias is of current research interest P|-p!3|]. Such possibihty requires the 
system to be out of equihbrium in order for the process to be consistent with the second law of 
thermodynamics. Several physical models have been proposed to obtain such motion. In all the 
models noise plays the central role. One of the most discussed models is the one in which an 
asymmetric periodic potential system is adiabatically rocked by applying constant forces F 
and —F at regular intervals of time. One obtains unidirectional motion because in such a system 
the current i{+F) ^ —j{—F). Eventhough the time averaged applied force over a period vanishes 
the averaged current < j >= 0.5[j{+F) + j{—F)] becomes finite in the presence of noise ( thermal 
fluctuations). Moreover, the average current < j > was found to peak at an intermediate noise 
strength (or temperature). In this model it has been further shown that by suitably choosing the 
asymmetric periodic potential one may obtain current reversal |TB[ as a function of temperature 



provided the rocking frequency is high. Similar results, however, can be obtained in the presence of 
a unbiased colored noise instead of the oscillating force. There are several other interesting models 
to obtain unidirectional motion including models where potential barriers themselves are allowed 
to fluctuate or models wherein symmetric potential system is driven by temporally asymmetric 



forces pHllI, etc. 

The result that thermal noise helps to obtain unidirectional current in a periodic system was 
quite important. But later on it was pointed out that obtaining mere current does not necessarily 



mean that the system does work efficiently |]T^. Doing work involves flow of current against load 
and hence one must, in the spirit of the model, obtain current up against a tilted (but otherwise 
periodic) potential system. Analysis shows, however, that the efficiency of an adiabatically rocked 



system (ratchet) monotonically decreases with temperature. Therefore though such a ratchet system 
helps extract a large amount of work at an intermediate temperature ( where the current peaks) 
the work is accomplished at a larger expense of input energy; thermal fluctuation does not facilitate 



efficient energy conversion in this model ratchet system. In a subsequent work [|T^] this deficiency 
was rectified but in a different model wherein the asymmetric potential oscillates in time, instead 
of its slope being changed (rocked) between +F and —F adiabatically. In both these models the 
friction coefficient was constant and uniform in space. The present work makes a detailed study 
of the rocked ratchet system with nonuniform friction coefficient which varies periodically in space 

IMlllIlHlll- 



In this work we take the friction coefficient to vary with the same periodicity as the potential 
but with a phase difference, 0. The phase difference 0, the amplitude A of variation of friction 
coefficient, the amplitude Fq of rocking, the load, etc. affect the functioning of the ratchet in an 
intricate and nontrivial manner. The two of the important results we obtain are: (1) The efficiency of 
the adiabatically rocked ratchet shows a peak as a function of temperature, though the peak (which 
may or may not exist in case of spatially asymmetric potentials) position does not coincide with the 
temperature at which the current peaks, and (2) the current could be made to reverse its direction as 
a function of noise strength and the amplitude Fq even at low frequencies of rocking. These attributes 
are solely related to the medium being inhomogeneous with space dependent friction. It is worth 
noting that the introduction of space dependent friction, though does not affect the equilibrium 
properties (such as the relative stability of the locally stable states), changes the dynamics of the 
system in a nontrivial fashion. Recently it has been shown that these systems exhibit noise induced 
stability, it shows stochastic resonance in washboard potentials without the application of external 



periodic input signal, ||23|j2^ and also unidirectional motion in periodic symmetric potential (non 



ratchet-like) systems. |T5|,|T],||,|9 



In the next section we describe our model and obtain an expression for current and efficiency 
in the quasi-static hmit. In Sec. Ill we present our results. 



II. EQUATION OF MOTION IN INHOMOGENEOUS SYSTEMS 

The nature of correct Fokker-Planck equation in the presence of space-dependent diffusion 
coefficient (inhomogeneous medium) was much debated earlier. Later on the correct expression was 
found from a microscopic treatment of system-bath coupling. The motion of an overdamped particle, 
in a potential V{q) and subject to a space dependent friction coefficient 7(g) and an external force 



field F{t) at temperature T is described by the Langevin equation [15,23-25,28 



where ^{t) is a randomly fluctuating Gaussian white noise with zero mean and correlation : 
< ^(t)^(t') >= 26(t — t'). Here < .. > denotes an ensemble average over the distribution of the 
fluctuating noise ^(t). The primes in Eq. (1) denote the derivative with respect to the space variable 
q. It should be noted that the above equation involves a multiplicative noise with an additional 
temperature dependent drift term. The additional term turns out to be essential in order for the 
system to approach the correct thermal equilibrium state. We take V{q) = Vo(g) + qL, where 
Vo(g + 2n7r) = Vo(g) = —Vsin{q), n being any natural integer. L is a constant force ( load) 
representing the slope of the washboard potential against which the work is done. Also, we take the 
friction coefficient 7(g) to be periodic : 

7(g) = 7o(l — Asin(g + 0)), where is the phase difference with respect to Vo(g). The equation of 
motion is equivalently given by the Fokker-Planck equation 

dPjq^t) d 1 r^ dP{q,t) p/.^^p/ fo\ 



This equation can be solved for the probabihty current j when F[t) = Fq = constant, and is given 



. ^ ksTjl - exp(-27r(Fo - L)/kBT)) 

' t exp( /,^+^- 7(^) exp ^°(-);(;o-.).) U 

In the presence of space dependent friction and the phase lag 7^ 0, vr and in the absence of 
load j(-Fo) 7^ ~j(~-^o) even for a spatially periodic symmetric potential. Thus when the system 
is subjected to an external ac field F{t) the unidirectional particle flow (or rectification of the 
current) takes place. The phase lag cf) brings in the intrinsic asymmetry in the dynamics of the 
system. When an externally applied ac force changes slowly enough (quasi-static or adiabatic limit) 
i.e., when the time scale of variation of F{t) is much larger compared to any other time scales 
involved in the system we can readily obtain an expression for the unidirectional current. For 
a field F{t) of a square wave amplitude Fq, an average current over the period of oscillation is 
given by, < j > = | [j(-^o) + j(~-^o)]- This particle current can even flow against the applied 
load L and thereby store energy in useful form. In the quasi-static limit following the method of 



stochastic energetics |3j] it can be shown ||17| , ^8[| that the input energy Em ( per unit time) and 



the work W ( per unit time) that the ratchet system extracts from the external noise are given by 
Ein = ^Fo[j{Fo) — j(— Fo)] and W = |L[j(Fo) + j(— Fq)] respectively. Thus the efficiency ( rj) of the 
system to transform the external fluctuation to useful work is given by 

^_ L[m)+j{-F,)] 
^ Fo[m)-j{-Fo)y 

Henceforth all our variables like < j >, E^n, W, Fq, T are made dimensionless. The amplitude 
of potential V is set to unity as all other energy scales are scaled with respect to V. We evaluate 
< j >,W, Ein, and r] numerically |^ using Eq. (3). 



III. RESULTS AND DISCUSSIONS 

First, we present our results for average (net) unidirectional current in a symmetric periodic 
potential induced by adiabatic rocking and in the absence of load. We emphasize here that to 
obtain these currents the system must be inhomogeneous. The phase lag (except for = 0, vr, for 
which unidirectional current is not possible) plays an important role in determining the direction 
and magnitude of < j >. For 27r > > vr, in the presence of external quasistatic force F{t) and in 
the absence of load L, we have a forward moving ratchet (current flowing in the positive direction) 
and for < </) < TT we have the opposite. For instance, if we examine the effect of friction coefficient 
close to the minimum of the potential two different situations are encountered depending on the 
value of (j). When > vr and Fq > the particle experiences lower friction near the barriers in the 
direction of acquired velocity. The situation is reverse when Fq < 0. From Eq. (3) it follows that in 
a static force Fq, j(-Fo) 7^ ~j (~-^o) , hence rectification of current occurs in the presence of external 
adiabatic drive. Moreover, it should be emphasized that the magnitude of current or mobility, in 
the static field Fq, depends sensitively on the potential and the frictional profile over the entire 
period. Depending on the system parameters the current or mobility can be much larger or smaller 
than the current or mobility of a particle moving in a homogeneous medium characterised by the 
space averaged frictional coefficient. In fact, in the intermediate values of temperature and Fq the 
mobility can be made much larger than their asymptotic limits. This leads to stochastic resonance 
in a washboard potential in the absence of ac signal [^,0 . 



In Fig. 1 the average unidirectional current < j > is presented as a function of (p and T. For 
this figure the value of amplitude of square-wave ac field Fq = 0.05, and the amplitude of frictional 
modulation A = 0.9. It can be seen that < j > changes sign as a function of and current exhibits 
either a minimum or a maximum as a function of temperature depending on the value of 0. At the 



two limits of temperature (0 and inf) the currents vanish. This is the case for the value of Fo less 
than the critical value Fc of Fq where the barrier to motion in either direction vanishes. In our case 
the critical value of Fq is equal to 1. This stochastic resonance- like phenomenon has been observed in 
rocked ratchet systems characterised by asymmetric periodic potentials |]3|j5|JT6|] . From the contours 
of the plot it is clear that as we move away from phase shift = tt positive (or negative) direction 
the temperature at which maxima (minima) occurs shifts to a larger value and the absolute value 
of the current at the peak decreases. However, the present symmetric potential situation does not 
lead to multiple current reversals as a function of temperature in the quasi-static limit of an external 
drive. 

The current as a function of T and Fq is shown in Fig. 2 for = l.Svr and A = 0.9. For 
smaller fields Fq compared to the critical field F^. the current exhibits a maximum as a function of 
temperature. As we increase Fq the temperature at which the peak occurs decreases. For fields 
larger than Fq the current, however, decreases monotonically with temperature because the barrier 
to motion disappears. In this high field region mobility of a particle decreases with increase of 



temperature However, as we increase the field the net current < j > monotonically changes 
and saturates to a finite value. This is in contrast < j > vanishes for large Fq in the ratchet subjected 
to a rocking force in the absence of space dependent friction. In Fig. 3 we have plotted < j > versus 
< Fq > and for fixed values of T = 0.5 and A = 0.9. It can again be seen clearly that the current 
monotonically varies and saturates to a value given by —^sin{(l)) independent of temperature. This 



result ||2^ follows from the analysis of Eq.(3). As expected current reversal can be seen as a function 
of for large Fq. It can also be verified that the current < j > increases monotonically as a function 
of the amplitude A of the friction coefficient and hence we do not present variation of < j >, etc., 
with A. 

We now discuss the efficiency of a symmetric periodic potential system in the presence 



of space dependent friction driven by an adiabatic periodic field. Tlie efficiency of such a system 
with uniform friction has been studied earher and it has been shown that temperature does not 



facihtate the efficiency {r]) of energy conversion in the system ^7^. To calculate rj we make use 
of Eqs. (3) and (4). In our analysis the load L is applied against the direction of net current (in 
the absence of load). In this situation particle current can flow against the applied load L less 
than some critical value Lc thereby storing energy in useful form. For L > one cannot talk 
meaningfully the concept of efficiency as the current flows in the direction of the load and hence 
no storage of useful energy takes place. In Fig. 4 we have plotted efficiency r], input energy Ein 
and work done W (scaled up by a factor 60 for convenience of comparison) as a function of T for 
the parameter values, Fq = 0.5, = l.Svr, A = 0.9, and the load L = 0.04. The figure shows that 
the efficiency exhibits a maximum as a function of temperature indicating that thermal fluctuation 
facilitates energy conversion. This in contrast to the case of uniform friction coefficient where rj 
decreases monotonically with the increase of temperature in the same adiabatic limit 0]. It is to 
be mentioned that the temperature corresponding to the maximum efficiency is not the same as the 
temperature at which the average current < j > becomes maximum in the absence of load. The 
temperature at which the extracted work maximizes is not the same as the temperature at which 
the efficiency becomes maximum for the same parameter values. The input energy increases with 
temperature monotonically and saturates at the high temperature limit, t], W, and Ein show similar 
qualitative behaviour for other parameter values. The above important observation of temperature 
facilitating the energy conversion is applicable for the spatially symmetric potential. In general in 
adiabatically rocked systems with frictional nonuniformities the increasing thermal noise need not 
increase the efficiency. The efficiency is sensitive to the qualitative nature of the periodic potential 
(and also to the nonuniformity of friction). For instance, asymmetric potential exhibits quite complex 
behaviour of rj and < j >. To illustrate this we take Vo{q) = — sin q — j sin 2q ( where /i lies between 



-1 and 1, and is the asymmetry parameter). With this potential we discuss three separate cases: 
Case A - system in a symmetric potential (/i = 0) in an inhomogeneous medium (A 7^ 0), Case B 
- system in an asymmetric potential (/x ^ 0) in an inhomogeneous medium (A 7^ 0), and Case C - 
system in an asymmetric potential [p, ^ 0) in a homogeneous medium (A = 0). 

In Fig. 5, we have presented results of 77 versus T for all the three cases described above. 
For this we have taken F — 0.5, L — 0.002, and = I.Stt. For case A A = 0.9, for case B = 0.08 
and A = 0.9, and for case C //. = 0.08. As discussed earlier for the case A temperature maximizes 
the efficiency. Case C, where the medium is homogeneous, efficiency monotonically decreases with 
temperature. These observations have been emphasized in earlier literature. The case B, where 
potential asymmetry and frictional inhomogeneity are present, the efficiency decreases monotonically 
in this parameter regime. In general whether the temperature facilitates the energy conversion in 
case B depends sensitively on the system parameters. In some limited parameter range the peaking 
behaviour is seen as in Fig. 6. The parameter values are indicated in the caption. The presence of 
asymmetry in the potential may or may not help in enhancing the efficiency of the system. This 
can be seen from Figs. 5 and 6. In all these cases the work done W and the input energy 
show similar qualitative features as shown in the inset of Fig. 1. Now, we discuss the variation of 
efficiency as a function of load. 

On general grounds it is expected that the efficiency too exhibits maximum as a function of 
load. It is obvious that the efficiency is zero when load is zero. At the critical value (beyond which 
current flows in the direction of the load) the value of current is zero and hence the efficiency vanishes 
again. In between these two extreme values of load the efficiency exhibits maximum. Beyond L — 
the current flows down the load and therefore the idea of efficiency becomes invalid. In Fig. 7, we 
have plotted r\ versus load for all the three cases for chosen values of parameters as mentioned in 
the figure. In all these cases current monotonically decreases as a function of load. The work done 

9 



against load W exhibits a maximum as a function of load. The load at which W shows maximum 
does not coincide with the load at which t] becomes maximum. The input energy Ein as a function 
of load varies non monotonically exhibiting a minimum. However, depending on the case under 
consideration the value of the load at which the minimum in the input energies observed may be 
larger than Lc above which efficiency is not defined. 

In Fig. 8, we have plotted the efficiency versus the amplitude of the adiabatic forcing Fq 
for all the three cases. It can be seen from the figure that for the system in an inhomogeneous 
medium, namely for cases A and B, t] exhibits a maxima and saturate to the same value in the large 
amplitude limit. In contrast, for the case C after exhibiting maximum rj goes to zero. This follows 
from the simple fact that in the large amplitude limit in the absence of frictional inhomogeneities 
the net unidirectional current tends to zero. The peculiar feature of saturation of efficiency in 
inhomogeneous media is related to the fact that the average current saturates to a constant value 
in the high amplitude limit as discussed earlier. This somewhat counter-intuitive result is typical to 
inhomogeneous media. Having discussed efficiency of energy conversion we now study the nature of 
net current < j > in the presence of spatially asymmetric potential to examine if current reversals 



take place in the adiabatic limit in the absence of load. It is known from the earlier literature |jT6 
that in an adiabatically rocked asymmetric potential ratchet system net current does not exhibit 
reversals as a function of T. In these systems current reversals are possible when the frequency 
of the applied ac field is large. We show here that in the presence of frictional inhomogeneities in 
addition to asymmetry in the potential one can observe current reversal as a function of thermal 
noise. In Fig. 9, we have plotted the magnitude of net current < j > versus T for all the three cases 
A, B, and C. The corresponding parameter values are mentioned in the caption of the figure. The 
cases A and C do not exhibit current reversal. This is a general result independent of parameter 
values. However, in case B current reversal is observed. To obtain current reversal both asymmetry 
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in potential and nonuniform friction coefficient are essential, that is, current reversals arise due to 
the combined effect of 4> and /i. Moreover, it should be noted that to observe current reversals the 
parameter range should be such that the net current in case A is in the opposite direction to that 
in case C. For the case B for which current reversal is observed, the plot of efficiency separates into 
two disjoint branches as the load should be reversed keeping the magnitude same when the current 
reversal takes place. In the presence of 7^ and A 7^ 0, where the current reversals are observed, 
the efficiency as a function of temperature is less than the maximum value of efficiency in either 
of the two cases A and C. That is, rjdj, ^ 0, A 7^ 0) < max[r){^i = 0, A 7^ 0),77(// 7^ 0, A = 0)]. To 
further analyze the nature of current reversals, in Fig. 10, we have plotted < j > as a function of 
(f) and T for fixed values of/x=l, A = 0.9, and Fq — 0.5. From the contour plots it is clear that 
as a function of the current reverses sign twice in the intermediate temperature range. Thus the 
current exhibits reentrant behaviour as a function of phase which is special to the case B. As we 
decrease the asymmetry in the potential the < j >= contour line shifts towards T = thereby 
enhancing the domain of current reversal to a lower value of temperature as a function of 0. As a 
function of temperature the current reversals occur in a definite range of phase which , in turn, 
depends on other material parameters. The qualitative behaviour of < j > remains unaltered for 
different Fq as long as Fq is less than the critical value. 

In Fig. 11, < j > is plotted as a function of T and Fq for = 1, = I.Stt, and A = 0.9. 
As opposed to the case of symmetric potential (Fig. 2), currents in the small temperature regime 
do exhibit maxima and then saturate to a constant value as noted earlier. As a function of T the 
current exhibits similar features as in the case A (Fig. 2). Figure 12, shows < j > as a function of 
Fq and for T = 0.1, /i = 1, and A = 0.9. As opposed to case A (see Fig. 3) < j > shows current 
reversal as a function of Fq in the range < < tt. However, in the asymptotic limit of Fq current 
saturates to a value — ■|sin(0) independent of the value of the asymmetry parameter ji. From the 
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contour plot it follows that as a function of phase cf) we observe the reentrant behaviour of current 
at high values of Fq. As we decrease ix the < j >= contour shifts towards smaller values of Fq 
thus making it possible to observe the double reversals at even smaller values of Fq. This reentrant 
behaviour as a function of (j) and the current reversal as a function of Fq is very specific to the case 
B alone (compare Fig. 3). 

Thus we conclude from our studies that the dynamics of a particle in an inhomogeneous 
medium is rich and complex. In the presence of adiabatic forcing and asymmetry in the potential 
current reversals can be observed as a function of T and Fq. And depending on the system param- 
eters thermal fluctuations facilitate the energy conversion. The above behaviour cannot be seen in 
the homogeneous medium in the same adiabatic limit. However, it is possible to observe these in 
homogeneous media in the presence of finite frequency ac drive (nonadiabatic regime). This seems 
to suggest that cf) may play the characteristic role of frequency in our model in the absence of nona- 
diabatic ac drive. This has been noted earlier in the context of observation of stochastic resonance 



phenomena in inhomogeneous media in the presence of static tilt alone p4[ . As a function of phase 
0, we observe reentrant behaviour for the current which arises because of interplay between asym- 
metry, inhomogeneity, thermal noise, and strength of the adiabatic forcing. Some of the phenomena 
can be understood at best at a qualitative level only. The effect of nonadiabatic forcing may be of 
further interest. Work on this line is under investigation. 
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FIGURE CAPTIONS 



Fig. 1. The net current < j > as a function of and T, for parameter values Fq — 0.5, 
A = 0.9. In the base plane contour of surface plot < j > (0, T) are given, dotted line indicates 
< j >— contour. 

Fig. 2. < J > as a function of temperature T and the amplitude of the rocking force Fq, 
for = I.Stt, a = 0.9. 

Fig. 3. < J > as a function of the amplitude of the rocking force Fq, and phase (p, for 
T — 0.5. In the base plane contour of surface plot is given. 

Fig. 4. Efficiency rj, and 1^ as a function of T for = I.Stt, Fq — 0.5, A = 0.9, and 
L = 0.4. W has been scaled up by a factor 60 to make it comparable with rj and Ei^- Y-axis is in 
dimensionless units. 

Fig. 5. Efficiency versus temperature for (i) case A (// = 0, and A = 0.9), (ii) case B 
{H = 0.08 , A = 0.9), and (iii) case C (/^ = 0.08 , A = 0.0), for fixed Fq = 0.5, = I.Stt, and 
L = 0.002. 

Fig. 6. Efficiency 77 versus temperature T, for different values of ji and for Fq = 0.5, A = 0.9, 
L — 0.1, and = I.Stt. 

Fig. 7. Efficiency versus load L for (i) case A (// = 0, and A — 0.9), (ii) case B (// = 1.0 , 
A = 0.9), and (iii) case C (/x = 1.0 , A = 0.0), for fixed Fq = 0.5, = I.Stt, and T = 0.1. 

Fig. 8. Efficiency versus Fq for (i) case A (/x = 0, and A = 0.9), (ii) case B [ji — 1.0 , 
A = 0.9), and (iii) case C (// = 1.0 , A = 0.0), for fixed L = 0.02, = I.Stt, and T = 0.1. 
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Fig. 9. Current < j > versus temperature T for (i) case A {/j, — 0, and A = 0.9), (ii) case B 
(// = 1.0 , A = 0.9), and (iii) case C (^i = 1.0 , A = 0.0), for fixed Fq = 0.5, = O.Stt, and L = 0.0. 

Fig. 10. The net current < j > as a function of 4> and T, for parameter values Fq — 0.5, 
A — 0.9, and /j, — 1.0. In the base plane contour of surface plot < j > {4>,T) are given, dotted line 
indicates < j >= contour. 

Fig. 11. < j > as a function of temperature T and the amphtude of the rocking force Fq, 
for (p = 1.3n, A = 0.9, and n = 1.0. 

Fig. 12. < j > as a function of the amplitude of the rocking force Fq, and phase (p, for 
T = 0.5, and /j, — 1.0. In the base plane contour of surface plot < j > {(f) ,T) are given, dotted hne 
indicates < j >= contour. 
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